Fidelity approach to Gaussian transitions 
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The fidelity approach to the Gaussian transitions in spin-one XX Z spin chains with three different 
values of Ising-like anisotropy A is analyzed by means of the density matrix renormalization group 
(DMRG) technique for systems of large sizes. We find that, despite the success in the cases of 
A = 2.59 and 1, the fidelity susceptibility fails to detect the Gaussian transition for A = 0.5. Thus 
our results demonstrate the limitation of the fidelity susceptibility in characterizing quantum phase 
transitions, which was proposed recently in general frameworks. 
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I. INTRODUCTION 



Due to latest progress in quantum information sci- 
ence people attempt to investigate the quantum 
phase transitions (QPTs) [2fl in many-body systems from 
the perspective of quantum information [3(. Recently, 
it was suggested that the ground-state fidelity, a ba- 
sic notion that emerged from quantum information sci- 
ence, should be able to signal the QPTs [J, Q. Be- 
cause the fidelity is a measure of distinguishability be- 
tween two states, one anticipates that the ground-state 
fidelity typically drops abruptly at the quantum criti- 
cal point, as a consequence of the dramatic changes in 
the structure of the ground states when systems un- 
dergo a QPT. Thus the fidelity should be able to fur- 
nish a signature of the quantum critical point. Another 
related transition indicator is provided by the singular- 
ity in the second derivative of the fidelity @, 0, HI (or 
the so-called "fidelity susceptibility" ((J). The fidelity 
approach has been tested for various man y- b ody sys- 

, [24, 



I. 



Since the fidelity is a purely Hilbert-space geometrical 
quantity, no a priori knowledge of the structure (order 
parameter, correlations driving the QPTs, topology, etc.) 
of the considered system is required for its use. There- 
fore, the fidelity analysis might be a potential universal 
tool for characterizing the QPTs. That is the reason 
why the investigation of the role of fidelity in detect- 
ing the topological phase transition, where no local order 
parameter can distinguish the ground states around the 
critical points, has started quite recently [29l. l30l l3l[ |32| . 

Besides the examination of its validity in restricted sit- 
uations, general understanding of the fidelity approach 
has come out, and the limitation of the fidelity suscep- 
tibility in detecting QPTs has been established [h| [2l[ . 
Conventionally, QPTs are characterized by the nonan- 
alyticities of the ground-state energy: first-order QPTs 



(lQPTs) are characterized by discontinuity in the first 
derivative of the ground-state energy, second-order QPTs 
(2QPTs) are characterized by discontinuity /singularity 
in the second derivative, and so on. It is pointed out that, 
while the fidelity susceptibility is indeed an effective tool 
in detecting the critical points of IQPTs and 2QPTs, it 
may fail to detect higher-order QPTs [2l|. Moreover, the 
authors in Ref. [19( show that, for critical (gapless) sys- 
tems of finite size L, the fidelity susceptibility S fulfills 
scaling relations 
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where d is the spatial dimension, z is the dynamic expo- 
nent, and Ay is the scaling dimension of the transition- 
driving term in the Hamiltonian. The scaling relation in 
Eq. ((T|) implies that, for critical points with Aq > 0, S 
is non-singular even in the thermodynamic limit. There- 
fore, the fidelity susceptibility is not suitable to single out 
these QPTs. 

Due to the lack of knowledge of exact ground state 
wave functions, the ground-state fidelity is usually diffi- 
cult to be calculated. For models that are not exactly 
solvable, most of researchers resort to numerical exact 
diagonalization on small systems. However, because of 
finite-size effects, the relevance of their results for the 
thermodynamic limit may not be guaranteed. There- 
fore, further investigations for systems of large sizes are 
clearly called for. By using the density matrix renormal- 
ization group (DMRG) technique [33[ , the fidelity suscep- 
tibility of the spin-one XXZ spin chains with an on-site 
anisotropic term for the case of the Ising-like anisotropy 
parameter A = 1 is evaluated for system sizes up to 
L = 160 [HI . The scaling relation in Eq. ([1]) has been 
confirmed numerically for both the Ising and the Gaus- 
sian transitions of this model. 

In the present article, we further pursue the fidelity ap- 
proach by analyzing the Gaussian transitions of the same 
model for various parameters of Ising-like anisotropy up 
to L = 400. Because the critical exponents of the Gaus- 
sian transitions change continuously along the critical 
lines, the QPTs of different order can be realized sim- 
ply by choosing different Ising-like anisotropy parameter 
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A. Here, we consider three cases of A = 2.59, 1, and 0.5, 
which correspond to 2QPT, third-order QPT (3QPT), 
and fifth-order QPT (5QPT), respectively (see Sees. M 
and IIIip . We find that, while the fidelity susceptibility 
can serve as a valid transition indicator for the Gaussian 
transitions in the cases of A = 2.59 and 1, it does not show 
any singularity around the critical point for A = 0.5 and 
fails to detect this higher-order QPT. Therefore, our re- 
sults provide a concrete illustration for the failure of the 
fidelity susceptibility in characterizing QPTs, and lend 
further numerical supports of the general proposals in 
Refs. [19] and [H[. 

This paper is organized as follows. The model Hamil- 
tonian and its low-energy effective theory, as well as the 
general scaling arguments, are described in Sec. HH Our 
DMRG calculations for the fidelity susceptibility are pre- 
sented in Sec. IIII1 We summarize and conclude our re- 
sults in Sec. IIVI The appendix contains an outline of the 
computation of the ground-state fidelity under DMRG 
algorithm. 



II. SPIN-ONE ANISOTROPIC MODEL AND ITS 
EFFECTIVE THEORY 

The S = 1 XX Z chains with uniaxial single-ion- type 
anisotropy is defined by the Hamiltonian: 
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where Sf (a = x, y, z) stand for the spin-one opera- 
tors at the j'-th lattice site. A and D parametrize the 
Ising-like and the uniaxial single-ion anisotropics, respec- 
tively. The full phase diagram consists of six different 
phases [H], [HI (see Refs. [H, H3, [H, H^j for recent nu- 
merical determinations). Here we focus our attention on 
the Gaussian transitions between the Haldane and the 
large-D phases. It has been found that the low-energy 
effective continuum theory for the Gaussian transitions 
can be described by the sine-Gordon model [37, 38] (here 
we follow the notations used in Ref. [38|) 
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n 2 + (d x $>) — j cos 



(3) 



where n and $ are the conjugate bosonic phase fields, 
and a is a short-distance cut-off of the order of the lat- 
tice spacing. The coefficient p, oc (D — D c ) in the vicinity 
of the critical point D c for a given A, and thus becomes 
zero at the transition point. The value of the Luttinger 
liquid parameter K varies continuously between 1/2 and 
2 along the critical line. Because all the scaling dimen- 
sions and the critical exponents are determined by a sin- 
gle parameter K, they also change continuously along 
the critical line. From the sine-Gordon theory [40j ], it is 
found that the critical exponent of the correlation length 
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and the scaling di mensio n Ay = K for the transition- 
driving term cos(V 4ttK^>). Thus the exponent in the 
scaling relation of the fidelity susceptibility in Eq. (fTJ 
becomes (d = 1 and z = 1 are assumed here) 



A< 
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(5) 



From this expression, one realizes immediately that, for 
critical points with K > 3/2, Aq is positive and then S 
becomes non-singular even in the thermodynamic limit. 
That is, the fidelity susceptibility fails to single out the 
Gaussian transitions with K > 3/2. 

According to the conventional classification, a 1QPT is 
characterized by a finite discontinuity in the first deriva- 
tive of the ground state energy. Similarly, a 2QPT is 
characterized by a finite discontinuity, or divergence, in 
the second derivative of the ground state energy, assum- 
ing the first derivative is continuous. From the effective 
theory in Eq. ([3]) and the scaling hypothesis, the first 
derivative of the ground-state energy density e(D) with 
respect to D at a fixed A gives [38] 
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where 



From Eq. 
tained as 



p = vAy = 
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the second derivative of e(D) can be ob- 
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Therefore, the second derivative of e(D) shows a diver- 
gence for < p < 1 (or < K < 1). In this case, the 
Gaussian transition corresponds to a 2QPT. Similarly, 
for 1 < p < 2 (or 1 < K < 4/3), the Gaussian transition 
can be classified as a 3QPT. Finally, as mentioned above, 
for the Gaussian transitions with K > 3/2, the fidelity 
susceptibility fails to detect them. In these cases, their 
values of p will be greater than three, and these transi- 
tions will be identified as the QPTs of higher than forth 
order. 



III. DMRG RESULTS 

In the present work, we consider three cases of A = 
2.59, 1, and 0.5, whose Luttinger liquid parameter has 
been found to be K = 0.85, 1.328, and 1.580, respec- 
tively [13, Hl| ■ From the discussions in the previous sec- 
tion, the QPTs of these three cases should correspond to 
2QPT, 3QPT, and 5QPT, respectively. Therefore, one 
expects that the fidelity susceptibility will not show sin- 
gularity for the case of A = 0.5 only. 
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In the following, our DMRG results are presented in 
order. The fidelity susceptibility for system in Eq. © of 
size L is calculated by [ljj, EH 



S(D) = lim 

8->V 



2[1 -T(D,D + S)} 
LS 2 



where the ground-state fidelity is given by [5j 
F{D,D + 5) = \(* (D)\V (D + 6))\ 



(9) 



(10) 



with \^q(D)) and l^ol-D + S)) being two normalized 
ground states corresponding to neighboring Hamiltonian 
parameters. Details on the computation of the ground- 
state fidelity under DMRG algorithm is explained in the 
appendix. In our calculations, S — 10~ 3 is used. Our re- 
sults are evaluated by means of the finite-system DMRG 
technique under open boundary conditions for system 
sizes up to L = 400, where 300 states per block are kept 
and five DMRG sweeps are performed for the trunca- 
tion error being 10 -9 at most. Finally, in order to select 
a specific ground state and stabilize our calculations of 
the fidelity around the Neel phase, a Zeeman term hiSf 
acting only upon the spin at site 1 is added in the Hamil- 
tonian of Eq. @ • Here we take the boundary magnetic 
field h x = -1. 



A. case of A = 2.59 

The findings of the fidelity susceptibility S(D) and the 
ground-state fidelity T{D, D + S) for A = 2.59 are shown 
in Fig.QJa). As shown in the inset, drops in the ground- 
state fidelity are observed, which signal precursors of the 
Gaussian transition in the model under consideration. 
Further evidences for indicating QPTs are provided by 
the results of the fidelity susceptibility S and the second 
derivative of the ground-state energy density d 2 e/dD 2 . 
As seen from Fig. [TJ the maximum values of both S and 
—d 2 e/dD 2 grow with increasing size, and thus indicate 
divergence in the L — > oo limit (see also Fig. [3] below). 
As discussed in Sec. [Til the divergent behaviors both in 
the fidelity susceptibility and the second derivative of the 
ground-state energy density indicate that the transition 
we found should be a 2QPT. 



According to the finite-size scaling theory [4l| , one has 
\D ma , x (L) — D c \ oc LT y l v , (11) 

where D c is the critical point in the thermodynamic limit 
and v is the critical exponent of the correlation length. 
Thus D c can be determined by an extrapolation proce- 
dure from the locations D max (L) of the local maxima in 
S and —d 2 e/dD 2 on a size-L system. The results for 
A = 2.59 is shown in Fig. H We find that D c ~ 2.293 
and v ~ 0.79 for the findings of S, while D c ~ 2.294 and 
v ~ 0.77 for the data of —d 2 e/dD 2 . Because of the rela- 
tion in Eq. the Luttinger liquid parameter K ~ 0.74 
(K ~ 0.70) for the data related to S (~d 2 e/dD 2 ). We 
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FIG. 1: (Color online) (a) Fidelity susceptibility <S, (b) 
the second derivative of the ground-state energy density 
—d 2 e/dD 2 as functions of D for various sizes L with A = 2.59. 
Inset in the top panel shows the fidelity T(D, D + S) as func- 
tions of D for the corresponding sizes. Here we take 8 — 10~ 3 . 



find that our results are consistent with the previous find- 
ings [H, where D c ~ 2.30 and K ~ 0.85. We note that 
the DMRG calculations in Refs. (3?| and are under 
the periodic boundary conditions, rather than the open 
boundary conditions used in the present work. The small 
discrepancy of our results from theirs may be due to the 
different boundary conditions employed. 

Moreover, at the critical point, where the only length 
scale is provided by the system size itself, the scaling 
formula in Eq. ([5]) implies 



d 2 e(D) 



L -(p-l)/u 



dD 2 

From Eqs. (|4]) and |J7}, one has 

-(p -l)/i/ = 2(1 -A") 



(12) 



(13) 



On the other hand, the fidelity susceptibility S is pro- 
posed to fulfill the scaling relation in Eq. (p}. To verify 
the predicted critical scaling behaviors in Eqs. (Q]) and 
(fl"2)) . the values <5> max and [— d 2 e/dD 2 ] max of the local 
maxima in S and —d 2 e/dD 2 for various sizes L are plot- 
ted in Figs.[3{a) and (b), respectively. It is found that our 
data do fulfill the predicted critical scaling behaviors with 
A Q ~ -1.28 and -(p - l)/u ~ 0.40. From Eqs. © and 
(P2J), our values of Aq and —{p—V)jv gives K ~ 0.86 and 
0.80, respectively. The Luttinger liquid parameter K de- 
termined by the present finite-size scaling agrees with the 
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FIG. 2: Finite-size scaling of D max of 5 (x) and -d 2 e/dD 2 
(A) versus L" 1 ''" for A = 2.59. The lines are least square fits, 
where v ~ 0.79 (y ~ 0.77) for those D ma x's corresponding to 
the local maxima in the curves of S {—d 2 e/dD 2 ). 
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FIG. 3: The log-log plot of (a) 5 max and (b) [~d 2 e/dD 2 ] m ^ 
for various sizes L with A = 2.59. The lines are least square 
straight line fits. 
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FIG. 4: (Color online) (a) Fidelity susceptibility S, (b) 
the second derivative of the ground-state energy density 
—d 2 e/dD 2 as functions of D for various sizes L with A = 1. 
Inset in the top panel shows the fidelity T(D, D + S) as func- 
tions of D for the corresponding sizes. Here we take 5 = 10 -3 . 



B. case of A = 1 

For the case of A = 1, a numerical support of the use 
of the fidelity in detecting QPTs has been provided in 
Ref. [23[, and the scaling relation of the fidelity suscep- 
tibility in Eq. JTJ) has also been verified there. Here, we 
focus our attention on the comparison between the scal- 
ing behaviors of the fidelity susceptibility and the second 
derivative of the ground-state energy density, which is 
shown in Fig. |U As observed in Ref. [23j], our present 
data in Fig. QJa) do show the developing drops in the 
ground-state fidelity and the divergent behaviors in S. 
Thus the fidelity susceptibility does signal precursors of 
the Gaussian transition in the model under consideration. 
However, as seen from Fig. 0|b), the second derivative of 
the ground-state energy density does not grow with in- 
creasing size in the same parameter regime. It implies 
that this QPT should be the one of higher than second 
order. To the author's knowledge, the present model at 
A = 1 may be the first concrete example in the literature 
with a QPT of higher than second order, which can be 
singled out unambiguously by using the fidelity suscepti- 
bility. 



previous findings [37| and those determined by the criti- 
cal exponent v coming from the scaling in Fig. [2] Thus 
the fact that a single parameter K controls all the criti- 
cal exponents for the Gaussian transition is confirmed by 
our numerical results. 



C. case of A = 0.5 

Although the fidelity susceptibility is possible to de- 
tect the QPTs of higher than second order, as explained 
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FIG. 5: (Color online) (a) Fidelity susceptibility <S, (b) 
the second derivative of the ground-state energy density 
—d 2 e/dD 2 as functions of D for various sizes L with A = 0.5. 
Inset in the top panel shows the fidelity T(D, D + S) as func- 
tions of D for the corresponding sizes. Here we take 8 = 10 -3 . 



in Ref. [2l| . it may not always work. For example, as 
discussed at the end of Sec. HH the fidelity susceptibility 
will fail to detect those Gaussian transitions with the Lut- 
tinger liquid parameter K > 3/2, which give Aq > and 
correspond to the QPTs of higher than forth order. In 
the following, we show that the present model at A = 0.5 
provides an example for this case. 

Our DMRG results of the fidelity susceptibility S(D), 
the ground-state fidelity T{D,D + S), and the second 
derivative of the ground-state energy density d 2 e/dD 2 
for A = 0.5 are shown in Fig. O While developing drops 
in the ground-state fidelity are indeed observed in the in- 
set of Fig. EJa), the fidelity susceptibility and the second 
derivative of the ground-state energy density show no di- 
vergent behavior up to L = 400. That is, this Gaussian 
transition must be a QPT of higher than second order, 
which the fidelity susceptibility fails to detect. Further 
support for such a conclusion is provided by the find- 
ings of <S max of the local maxima in S for various sizes 
L, which is plotted in Fig. [5J It is found that maximal 
values of iS do tend to saturate, rather than diverge, in 
the L — > oo limit. Thus this Gaussian transition does 
provide a counterexample of the usefulness of the fidelity 
susceptibility as an effective tool in detecting QPTs. 

One may wonder if the non-singular behavior in the 
fidelity susceptibility shown in Figs. [5] and [5] is due to the 
possible finite-size effects, even though sizes up to L = 
400 have been considered. By contrast, the entanglement 
entropy is evaluated below. It has been found that the 
divergent character of the entanglement entropy can also 
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FIG. 6: The log-log plot of <S> max for various sizes L with 
A = 0.5. The dashed curve is the least square fit by using the 
formula <S max = <Soo — aL~ b with <Soo = 0.073, a = 0.23, and 
b = 0.75. 



show the existence of the QPTs |42j, |43j, |4J, |45j, |4jJ |47|, |48( . 
Within the present DMRG calculations, the singular be- 
havior of the entanglement entropy is indeed observed in 
the given parameter regime (see Fig. [7] below). Thus the 
desired Gaussian transition does appear there. This indi- 
cates that the system sizes we employed should be large 
enough and finite-size effects should play no important 
role in our calculations. 

Here we consider the entanglement entropy, or the von 
Neumann entropy of the reduced density matrix p R {D) 
of the right-hand block of L/2 contiguous spins 



£(D) = -Tr[ PR (D)\og 2 p R (D)} 



(14) 



Our DMRG results are presented in Fig. [7] It is found 
from Fig. [7ja) that peaks do grow as size L increases, 
which show the existence of the expected Gaussian tran- 
sition. Further evidence is provided by the findings of 
£max of the local maxima in £ for various sizes L, which is 
plotted in Fig.UJb). It indicates that £ max oc l°g2(^) an d 
thus £ m ax will diverge in the L — > oo limit. It is known 
that the slope in Fig. [7Jb) gives the value of c/6 un- 
der open boundary conditions [USB S H S3, 13 , 
where c is the central charge of the conformal field theory. 
Our data give c ~ 0.98, which agrees with the theoretical 
value c = 1 for the Gaussian transition [37j. Further- 



more, according to Eq. (jTTJ) , the corresponding critical 
points D c can be deduced from the locations D max (L) of 
the local maxima in £{D) on a size-L system, as shown 
in Fig. Etc). It is found that D c ~ 0.63 and v ~ 1.51. 
The value of v implies the Luttinger liquid parameter 
K ~ 1.34 due to the relation in Eq. g]). In Ref. 
D c ~ 0.65 and K ~ 1.580 are reported. Thus our results 
are again consistent with their findings. As mentioned 
before, the discrepancy of our results from theirs may be 
due to the different boundary conditions employed. In 
short, the missing QPT in the measurements of the fi- 
delity susceptibility is able to be captured by using the 
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FIG. 7: (Color online) (a) Entanglement entropy £ as func- 
tions of D for various sizes L. (b) The semi-log plot of £ max 
for various sizes L. The line is least square straight line fit. 
(c) Finite-size scaling of D max of £ versus L~ 1/V . The line is 
least square fit with v ~ 1.51. Here A = 0.5. 



tool of the entanglement entropy. 



Moreover, the critical point and some of the correspond- 
ing critical exponents are determined through a proper 
finite-size scaling analysis, and these values agree with 
the findings reported in the literature. Although it is 
successful in characterizing the QPTs for the cases of 
A = 2.59 and 1, the fidelity susceptibility is not able 
to detect the Gaussian transitions for A = 0.5. This 
observation can be explained by the general theories in 
Refs. pjj and [2l[. Thus our work provides an instance 
for the limitation of the use of the fidelity susceptibility 
for revealing QPTs. 

Finally, we remind that developing drops in the 
ground-state fidelity are observed around the critical 
point even for the case of A = 0.5 [see the inset of 
Fig. [5(a)]. This seems to indicate that the singularity 
of the fidelity may still be a valid transition indicator 
even for higher-order QPTs. The failure of its second 
derivatives (say, the fidelity susceptibility) may just im- 
ply that higher derivatives of the fidelity are necessary 
to signal some higher-order QPTs. Indeed, such a con- 
clusion can be deduced from the general arguments in 
Refs. [3 and ft was shown in Ref. that the 

fidelity susceptibility S (or the second derivative of the 
fidelity per site) can be associated with the first derivative 
of the density matrix po of the ground state [see Eq. (6) 
therein]. Thus S can not detect the higher-order QPTs 
resulting from the nonanalyticities in second (and even 
higher-order) derivatives of po- That is, these higher- 
order QPTs can be recognized only by the singular be- 
havior in higher-order derivatives of the fidelity. Besides, 
the singular part of S has been shown in Ref. to 
behave as S ~ \D — D c \' /Aq . Therefore, the third deriva- 
tive of the fidelity should behave as \D ~ D c \" Aq ~ 1 . This 
implies that, although the second derivative of the fi- 
delity (or the fidelity susceptibility S) does not diverge 
for the Gaussian transitions with 3/2 < K < 5/3 [thus 
< vAq < 1, see Eqs. (0} and ©], the third derivative 
does become singular at these critical points. That is, one 
does need the third derivative of the fidelity to single out 
these higher-order QPTs. The more complete and deeper 
understanding of the role played by the fidelity and its 
derivatives in the studying of the higher-order QPTs is 
an interesting issue which calls for future investigations. 



IV. CONCLUSIONS 



In summary, we present numerical analysis of the fi- 
delity approach to the Gaussian transitions in the spin- 
one XX Z spin chains in Eq. @ with three different val- 
ues of Ising-like anisotropy A. It is found that the fidelity 
susceptibility can detect the Gaussian transitions for the 
cases of A = 2.59 and 1. As in the case of A = 1 in- 
vestigated in Ref. [23| , the scaling relation of the fidelity 
susceptibility proposed in Ref. [191 ] is again verified by 
the present DMRG calculations for the case of A = 2.59. 
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FIG. 8: The configuration of wave functions in the superblock. 
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FIG. 9: The basis transformation between a block and an 
enlarged block. 



APPENDIX: COMPUTATION OF FIDELITY 
UNDER DMRG ALGORITHM 

The DMRG method had been used to calculate the 
overlap of two different ground states for a decade. For 
example, people had applied this technique to evaluate 
the exponent of the orthogonality catastrophe for the 
problem of a single impurity in a one-dimensional Lut- 
tinger liquid [49[ . Recently, this method is used to inves- 
tigate the scaling properties of fidelity susceptibility for 
the model in Eq. @ with A = 1 23]. For completeness, 
the algorithm within the DMRG framework of computa- 
tion of the ground-state fidelity (i.e., the overlap between 

I 



Remind that (a'\a) and (f3'\f3) are nothing but the matrix 
elements of the s x s identity matrix 1 in the subspaces 
of the local Hilbert spaces at sites L/2 and L/2 + 1, re- 
spectively. Besides, let's define 

[i|_iJms,ms = (ms\ms) , (A. 4) 

2 L 

[if _i]m E ,m B = (m E \m E ) , (A. 5) 

which serve as the matrix elements of the m x m "iden- 
tity matrices" [If ] and [if ,] of the system and the 

2 1 2 1 

environment blocks of length ~ — 1, respectively. There- 
fore, the expression of the overlap in Eq. (|A.3[) can be 
rewritten as the following compact from 

($|*) = (($||[i]||*)) (A.6) 

with 

[i] = [i£_i]®i®i®[lf_i]> ( A - 7 ) 

2 1 2 1 

where 11$)) and [["J/)) represent the column vectors of 
dimension s 2 m 2 with the components <f>ms,ct' ,m E an d 



two ground states corresponding to two different system 
parameters) is presented in this appendix. Similar discus- 
sions in terms of the formalism of matrix product states 
can be found at the end of Sec. 2.1 in Ref. (50j . 

Let |$) and \*Sf) be the normalized wave functions ob- 
tained from two different DMRG calculations with dis- 
tinct Hamiltonian parameters for systems of the same 
size. Assume that the same number of states m is kept 
in both DMRG calculations. Moreover, both wave func- 
tions are expressed in the same scheme, where the cor- 
responding system (and thus the corresponding environ- 
ment) blocks are of the same length. Therefore, these 
two wave functions can be written as (see Fig. [8]) 

|*)= Y *ms,a,l3,mB\ms)\a)\0)\mis) , (A.l) 



|*)= Y ®rn s ,a',0',m E \rhs)\a')\P')\fh E ) . (A.2) 

iris ,ct' ,0' ,fh,E 

Here \a) da')) and \0) (|/3')) denote the states in the lo- 
cal Hilbert space of dimension s at sites L/2 and L/2+1, 
respectively. |m<j) and \fhs) represent the m renormal- 
ized states of the system blocks, while \m E ) and \m E ) are 
those m renormalized states of the environment blocks. 
Thus the overlap between these two ground states be- 
comes 



(A.3) 

I 

^ms,a,ftmB) respectively. Now the remaining work of 
evaluation of fidelity reduces to find a general strategy of 
calculating the two m x m matrices [if,] and [if ]. 

2 1 2 1 

Note that the renormalized states of the environment 
blocks of length I can be expressed as (see Fig. |9|) 

\m E ) = J2 [°Fk a ,m' B ), mB \a)\m' E ) , (A.8) 

a,m' E 

\fh E ) = Y, [OFh,fn' E ),fn e mrh' E ) , (A.9) 

f3,m' E 

where [Of ]( a ,m' B ),m E and [df]{p,m' E ),m B denote the ma- 
trix elements of an sm x m transformation matrices [Of] 
and [Of] for two different DMRG calculations with dis- 
tinct system parameters. They are formed by m eigenvec- 
tors with the largest eigenvalues of the reduced density 
matrix of the environment blocks of length I. By using 
Eqs. (|A.8|) and (|A.9|> . one obtains 



($|*)= Y ^s, a ',f3',rn E (ms\ms}(a'\a}(P'\/3}(m E \m E }^ mStat ^ 

rhs ,ct/ ,m£;;ms ,a,/3,m£; 



8 



(m B \m E )= ([A £ ](/3,m' E )^ E )^/3|a)(m' i5 |m' B )[Of] (ct , m ^ ) , mE . (A.10) 



Following the similar definition in Eq. (|A.5|) . the above 
expression can be casted into a compact form: 

[if] = [of] t (i®[i^ 1 ])[of], (A.ti) 

which gives the recursive relations of the m x m "identity 
matrices" [if] of the environment block of length I. For 
instance, in the cases of / = 2 and 3, one has 

[if] = [of]t(i®i) [of], 
[if] = [o 3 B ]t(i®[if]) [oi]. 



Similar recursive relations for the system blocks reads 

[if] = [of] + ([if-i]®i)[of]- (A.12) 

If we iterate all the way from both ends of the sys- 
tems by employing the recursive relations in Eqs. (|A. 1 1[) 
and (|A.12|) , then the matrices [1?,] and [if J can be 

2 1 2 1 

constructed eventually. Substituting these results into 
Eq. (|A.6p . the value of the ground-state fidelity will be 
obtained. 
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